secutive radar profiles are offset slightly 
(by amounts much less than A r) with re- 
spect to each other according to the rel- 
ative change in geodetic altitude. As a 
consequence, in a case in which the sur- 
face area under examination is homoge- 
nous (e.g., an ocean surface), a se- 
quence of consecutive radar profiles of 
the surface in that area containssamples 
of the surface response with range reso- 
lution (Ap) much finer than the range- 
bin increment (Ap «A r). 

Once the high-resolution surface re- 
sponse has thus become available, the 
profile of surface clutter can be accu- 
rately estimated by use of a conventional 
maximum-correlation scheme: A trans- 


lated and scaled version of the high-res- 
olution surface response is fitted to the 
observed low-resolution profile. The 
translation and scaling factors that opti- 
mize the fit in a maximum-correlation 
sense represent (1) the true position of 
the surface relative to the sampled sur- 
face peak and (2) the magnitude of the 
surface backscatter. 

The performance of this algorithm has 
been tested on CloudSat data acquired 
over an ocean surface. A preliminary 
analysis of the test data showed a surface- 
clutter-rejection ratio over flat surfaces of 
>10 dB and a reduction of the contami- 
nated altitude over ocean from about 1 
km to about 0.5 km (over the ocean). 


The algorithm has been embedded in 
CloudSat LIB processing as of Release 04 
(July 2007), and the estimated flat sur- 
face clutter is removed in L2B-GEOPROF 
product from the observed profile of re- 
flectivity (see CloudSat product docu- 
mentation for detailsand performance at 
http:/ / www.cloudsat.cira.colostate.edu/ 
dataSpecs.php?prodid=l). 

This work was done by Simone Tanelli, 
Kyung Pak, Stephen Durden, and Eastwood Im 
of Caltech for N A SA'sJet Propulsion Laboratory. 

The software usd in this innovation is 
available for commercial licensing. Please 
contactKarina Edmonds of the California In- 
stitute of Technology at (626) 395-2322. 
Refer to NPO-44873. 


MODIS Atmospheric Data Handler 


Stennis Spa ceCen ter, M ississippi 

A number of science data sets are derived 
from the observations of the Moderate 
Resolution Imaging Spectroradiometer 
(MODIS) instrument onboard NASA's 
Terra and Aqua satellites. These data 
typically contain information on re- 
trieval techniques, quality-control flags, 
and geo-referencing information. 
These datasets, distributed in H DF (FI i- 
erachical Data Format), must be fur- 
ther processed to extract relevant infor- 
mation for weather analysis studies and 
numerical models input. The MODIS- 
Atmosphere Data H an d I er software 
converts the FI DF data to ASCII format, 
and outputs: (1) atmospheric profiles 


of temperature and dew point and (2) 
total precipitable water. Quality-control 
data are also considered in the export 
procedure. 

The package currently consists of pro- 
grams to process the MOD05 and 
MOD07 data products from MODIS. 
The software is written using the C pro- 
gramming language and contains Make- 
files for easier compilation and installa- 
tion. The MODIS-ADFI software helps 
ease the overhead involved in data pro- 
cessing so that the numerical modelers 
may concentrate on their science and 
modeling tasks rather than manipulat- 
ing data for their models. 


This program was written by Valentine 
Anantharaj and Patrick Fitzpatrick of the 
Northern Gulf Institute at M ississippi State 
University for Stennis Space Cen ter. 

Inquiries concerning rights for its commer- 
cial use should be addressed to: 

M ississippi State University 
Northern Gulf Institute 
BLDG 1103, Room 233 
Stennis Space Center, MS 39529 
PhoneNo.: (228) 688-1157 
Fax: (228) 688-7100 
E-mail: val@gri.msstate.edu 
Refer to SSC-00267, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 


6»M ultibeam Altimeter N avigation U paate U sing Faceted 
Shape Model 

T he model is applicable to a body having almost any complex shape. 

NASA's Jet Propulsion Laboratory, Pasadena, California 


A method of incorporating informa- 
tion, acquired by a multibeam laser or 
radar altimeter system, pertaining to the 
distance and direction between the sys- 
tem and a nearby target body, into an es- 
timate of the state of a vehicle upon 
which the system is mounted, involves 
the use of a faceted model to represent 
the shape of the target body. I n the orig- 
inal intended application, the vehicle 
would be a spacecraft and the target 
body would be an asteroid, comet, or 
similar body that the spacecraft was re- 
quired to approach. The method could 


also be used in navigating aircraft at low 
altitudes over terrain that is rough 
and/or occupied by objects of signifi- 
cant structure. 

Fundamentally, what one seeks to 
measure isthe distance from the vehicle 
to the target body. The present method 
is the product of a generalization of a 
prior method of altimetry, in which the 
target body has a simple shape repre- 
sented by a spherical or ellipsoidal 
model. In principle, the estimate of dis- 
tance or altitude obtained by use of a 
multibeam altimeter can be more robust 


than that obtained by use of a single- 
beam altimeter, but if the surface of the 
target body has a complex and/ or irreg- 
ular shape, then it becomes more diffi- 
cult to define the distance and compute 
the distance from readings of a multi- 
beam altimeter. 

The faceted shape model of the pres- 
ent method facilitates the definition and 
computation of distance to a target ob- 
ject having almost any shape, no matter 
how irregular and complex. The use of 
faceted shape models to represent com- 
plex three-dimensional objects is com- 
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RANGE-M EASUREM ENT TABLE 


# 

Confidence 

M 

a 

N T = 

[N x , N y , N z ] 

K 

d 7 ‘ = 

ld x , d y , d z ] 

c T = 

[C x , c y , c z ] 

1 








2 








• 









The Laser or Radar Altimeter measures the magnitude of the vector (r) from the origin of each beam 
to the intersection of the beam with a target body represented by a multifaceted shape model. There 
are multiple beams, but only one is shown here for simplicity. The measurement data and associated 
model data are listed in the range-measurement table. 


mon in the computer-graphics literature 
and in the movie and video-game indus- 
tries. I n this method, the distance to be 
measured isdefined as the length of the 
vector ( p) from the center of mass of the 
multifaceted shape model to the center 
of mass of the vehicle, as depicted in the 
upper part of the figure. 

The state-update information derived 
from the most recent set of multibeam- 


altimeter measurements is listed system- 
atically in a range-measurement table 
(RMT), depicted in the lower part of the 
figure, in which the planar facets of the 
shape model are represented in Hesse's 
normal form. Each row of the table con- 
tains the data from one of the altimeter 
beams. The first column contains the 
rowindex (/), which isthecardinal num- 
ber of the affected beam. The second 


column contains a number, between 0 
and 1, representing the degree of confi- 
dence in the measurements. At the pres- 
ent state of development of the method, 
the confidence is taken to be either 0 
(signifying complete rejection) or 1 
(representing complete acceptance) of 
the data in the row. The third column 
contains the scalar range measurement 
|r| of the /th beam; the fourth column 
contains the standard deviation (a) of 
the range measurement. 

The fifth column contains the Carte- 
sian components [ A/ x , N y , A/ z ] of the 
transpose of the unit vector (N r ) normal 
to the model facet containing the inter- 
section of the /th laser beam with the 
surface of the target object. Typically, 
this intersection point is not known ex- 
actly and must be estimated, on the basis 
of tire current state estimate, by a previ- 
ously developed method that lies be- 
yond the scope of this article. The sixth 
column contains the facet constant, k 
(the perpendicular distance from the 
center of mass of the target body to the 
affected facet). The seventh column 
contains the Cartesian components [cf x , 
d y , d z ] of the unit vector along the /th 
laser beam. The seventh column con- 
tains the Cartesian components [ c x , c y , c z ] 
of the position vector from the center of 
mass of the vehicle to the origin of the 
/th laser beam. 

The entries in the RMT are mapped 
into a measurement equation for use by 
a Kalman filter that incorporates altime- 
try information into the final estimate of 
the state of a spacecraft or other vehicle 
maneuvering in the vicinity of a target 
body. The relative position vector, p, is 
part of the state vector that is updated by 
use of the Kalman filter. 

T his work was done by David S. Bayard, 
Paul Brugarolas, and Steve Broschart 
of Caltech for NASA's Jet Propulsion Labo- 
ratory. For more information, contact iaof- 
fice@jpl.nasa.gov. NPO-44428 
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